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ABSTRACT 


This paper deals with the simulation of dynamic model of permanent magnet 
synchronous generator (PMSG) in D-Q axes of the rotor rotating reference frame. The 
iron core losses and stray load losses of the machine are taken into account. The iron 
core losses are represented by iron core resistance connected in parallel with 
magnetizing inductance and then reflected into the stator side as a voltage drop to 
prevent increasing the number of differential equations in the model. The modified 
equivalent circuit can deal with all machine parameters without losing the accuracy of 
generator performance calculations. The modified equivalent circuit can be used as an 
efficient tool for analysis, design, and vector control algorithm of this type of generator, 
especially in renewable energy utilization. The model is executed by Matlab Simulink 
and very good results are obtained and compared with the results of the experimental 
model to display the validity and accuracy of the proposed dynamic model. 
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1. INTRODUCTION 


The main function of an electrical generator in renewable energy systems is to convert the 
rotational torque from wind or hydropower turbine as the prime mover into electrical power to 
feed the isolated load or the grid. Different kinds of generators are being used with hydro and 
wind turbines. The common type of AC generators that are candidates in micro-hydro and wind 
turbine systems are squirrel-cage induction generators and permanent magnet synchronous 
generators (PMSG). 
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Due to its advantages, recently, the PMSG is used in wind power systems. Its usefulness 
over the other kinds of generators are : 


e Due to the absence of brushes, rotor windings, and slip-rings, this type of generator has 
high reliability. 


e There is no additional excitation power supply due to the permanent magnet. 
e There is no rotor copper loss due to no rotor winding. 
e Lower maintenance, higher capacity, and higher massive/weight ratio. 


e This type of generator can be made with a large amount of magnetic poles and the 
turbine rotor coupled directly to the generator shaft. 


e Due to the low generator speed, there is no need for a gearbox. 


Because of the many benefits of this type of generator in renewable energy conversion, all 
researchers [1-23] mentioned in this work have taken up this topic to analyze and represent this 
type of generator employing a dynamic model. But all those who studied the mathematical 
model of the generator type (PMSG), used an incomplete dynamical model, as the iron core 
losses, stray load loss, and magnetizing saturation effect are not taken in the dynamic model. 
Also, there are no results in their analysis indicating the validity of the mathematical model. 
These types of machine losses cannot be omitted when performance and efficiency 
determination is required. Also, in all those dynamic models the excitation capacitance, which 
is linked in series with each stator phase of the generator is not considered. The excitation 
capacitance can be linked in series with machine terminals to obtain maximum output power 
from the generator at synchronous speed and resonance conditions. 


In the new and comprehensive dynamic model used in this research, the effect of machine 
stray load losses and iron core losses were included. The magnetizing saturation effect is also 
considered in the model due to this type of generator is operating in a non-linear saturation 
region at full load. The machine equivalent circuit is modified for comprehensive analysis 
without increasing the number of system differential equations or losing the accuracy of the 
performance calculations by converting the iron core resistance into a voltage drop on the stator 
circuit. The capacitor was also inserted in series with the machine terminals to get the higher 
output power from the generator at synchronous speed and resonance conditions. 


This work 1s structured in four sections; the first section 1s dealt with the introduction, the 
second section is dealt with the detailed dynamic model of the PMSG. The third section deals 
with the results and discussion while the last section deals with the conclusion. 


2. DYNAMIC MODELING OF PMSG 


The permanent magnet synchronous generator (PMSG) in figure (1) can generate electrical 
energy if its rotor rotating at a suitable speed by wind or hydro turbine as an external prime 
mover. The generator terminal voltage 1s build-up due to the permanent magnet fixed on the 
rotor with a suitable number of magnetic poles. Capacitors are linked in series with each phase 
to obtain higher output power from the generator at resonance condition. 
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3-Phase load 





Figure 1 PMSG driven by a prime-mover 


Figure (2) exhibits the per-phase equivalent circuit of the PMSG in D-Q axes rotating frame 
with inductive load and series capacitance, taking the effect of machine iron core losses and 
stray losses. The iron core losses and stray losses are represented by equivalent iron core 
resistance (R;,) connected in parallel with the magnetizing inductance, and stray load loss 
resistance (R,,,~) connected in series with the stator circuit. 
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Figure 2 D-Q axes equivalent circuit of PMSG 
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Figure (3) shows the modified equivalent circuit of the PMSG. In this circuit, the parallel 
branch of the iron core resistance (R;,) 1s replaced by series resistance and then reflected as a 
voltage drop in the stator circuit [24]. The magnetizing inductance in D-Q axes (Lyg and Lygq) 


are modified. 
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Figure 3 Modified equivalent circuit of PMSG in D-Q axes rotor rotating frame 


The stator series stray load resistance (R,,<p), the modified iron core resistance 
(Rssica ANd Rssicg) and the modified magnetizing inductance (Lyq and L,,,) in the D-Q axes 


model can be obtained as [25, 26]: 


Rssse — (W,. = Lge) * R,/[Ré i (W,, = Lg)" | (1) 
Rssica -_ Ric * Wr * Lua/(Rf, + (w,, ‘i Lua)‘] (2) 
Rssicq — Ric * Wy * Luq/|Ré: + (w,, ii Lug)" | (3) 
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Lima = Lua - Ri, [Ric P (w, - Lua)‘ (4) 
img = LMq * Ri, [Ric + (Wy * lag)? | (5) 
Where R, : the stator phase resistance. 
w, : rotor speed in electrical radian/second. 
Lua : D-axis magnetizing inductance/phase. 
Lyq : Q-axis magnetizing inductance/ phase. 
R;-: ron core resistance obtained from the no-load test as: 
Ric=Vo/l Pai - I *(RstRssse)] 
V, & I,: no-load phase voltage and current when the machine is driven as a motor. 
P,,,: no-load power per phase when the machine is driven as a motor. 


The D-Q axes series capacitance (C,gand C,,) in each phase must satisfy the following 
condition at resonance and synchronous speed to obtain maximum generator output power: 


Coq = 1/[we * (Lua + Les + Le)] (6) 

Coq = 1 /[W? * (Lug + Les + Le)| (7) 
Where Ly» : the load inductance per phase. In case of round or cylindrical rotor machine: 
Lua = Luq = Lu (8) 

Coa = Coq = C, (9) 


The voltage in D-Q axes of the stator and rotor of PMSG with series capacitance and 
inductive load derived from the modified circuit of figure (3) as: 


; / di, 
Vsq = —Vca — Rs * lsq — Rssse * lsq — Rssica * (i i” isa) + Wy * Lg * lsg — La ‘i —s (10) 
F 7 . F rf dis 

Vsq = —Veq — R, * lsg — Resse * sg — Rssica * log — Wy * La * log 1 Wy Pr = Lg ss os (11) 
Vrq = Wr Pr — lsq * Lma (12) 
Urq = —Wy * Ling * Isq (13) 
Vca = J Ge + Ww, * Veq) dt + Vedo (14) 

ls 

Veg = J (+ wy * Vea) dt + Vego (15) 
Psa = —Les * lsq + (i, — Leg) * Ling = —ha * tsq + db, (16) 
Pog = —(Les + Ling) *isq = hg * isq (17) 


Where Vsq 5 Vsq ,lsq ANd isg , Stator phase voltages and currents in D-Q axes respectively. 


, is the permanent magnetic flux per phase referred to the stator circuit? 


ir is the field current generated in the rotor due to the permanent magnet referred to the 
stator circuit? 


ip = br / Lima (18) 
La = (Les + Lmq) is the stator d-axis self-inductance. 

Lg = (Les + Lm) is the stator q-axis self-inductance. 

For round rotor machine Lg = Lg = Ls 


The stator (r.m.s) phase voltage can be obtained as: 


Von = |v2q + v2, |?/V2 (19) 
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In the case of resistive-inductive load the components of phase voltage (Vsq and Vsq) can 


be calculated as: 
dicg 








Vsq = Re * lq + Le * dt W, * Le * sq (20) 
: dl : 
Vsq = Re * isg + Le * = + wy * Lp * igg (21) 


Where Rp and Ly , load resistance and inductance per phase. 


The stator phase currents (i,q & i,q) and their derivatives in first-order differential equations 
can be obtained from equations (10) and (11) as: 


as " la : 3 
—s = (1/Lq) . |-Vsa — Vea — (R, a Rssse) *isq — Rssica * (i, — isa) + W, * Lg * Le (22) 


dis c 
ce (1/Lq) * |-Vsq — Veq — (Rs + Resse + Resicq) * isq — Wr * La * isa + Wr * br] (23) 


The electromagnetic torque generated by PMSG can be obtained as: 


Te = (=) (Py) * [Br * isqg — (La — Lg) * isa * icq] (24) 
In the case of a round rotor or cylindrical machine, the torque becomes: 
Te = (=) (P,) : by * Usq (25) 


Where F, magnetic pole pairs. 


The generator mechanical speed (w,,,) in radian/second can be obtained as [1]: 

dw 1 

Te = 7 Tm — Te — Kp * Wm) (26) 
Where J total moment of inertia in Kg.m”. 

T,, prime mover torque in N.m. 

Ky friction coefficient (N.m/(rad./sec. )). 


The rotor speed in electrical radian/second can be obtained as: 


W, = Py * Wm (27) 
The generator input power can be obtained as: 
P= Te Wy, (28) 
The generator output power delivered to the load can be calculated as: 
Pout = Pin — 3 * Loh * (Rs + Rssse) (29) 
Where [,,, load current per phase and can be obtained as: 
1 

Inn = [i2g + i2|?/V2 (30) 
Alternatively, the load active and reactive power can be obtained as: 

6) : ; 
Pp = (=) : (Vsa * lsq + Usq * isa) (31) 

3 5 ; 
Qe = (=) is Ce * log — Vsq * isq) (32) 
The load power factor can be calculated as : 

1 

P.F.= Pp/[Pz + Qe]? (33) 
The generator efficiency can be calculated as: 


n% = 100 (34) 


fa ee 
(Pint+Prot) 


Where P,o¢ 18 the total rotational (mechanical) losses for generator and prime mover. 
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The PMSG usually has higher efficiency than the other types of generators due to the fixed 
magnet for the generation of rotor flux. 


3. DYNAMIC SIMULATION RESULTS 


The dynamic simulation model of PMSG 1s developed in Matlab Simulink. 


The PMSG is driven experimentally at synchronous speed by a DC motor as a prime mover 
and the capacitors are connected in series with generator terminals to obtain maximum output 
power from the generator at resonance conditions. The parameters of the PMSG are given as : 


PMSG Parameters: 3-phase, round rotor, 230 V per phase; 2 kW with the following 


parameters : 
R,, = 3000; R, =30; Ry552 = Rp =200;L, =0.1 mH ; 
Lp, = 10 mH ; d, = 0.11 wb per phase ; P, =2 ; Ky = 0.003 N.m/(rad./sec. ) 
Jg = 9.03 kg .m* __. The dynamic magnetizing inductance which is found from the no- 


load test in terms of magnetizing current (J,,), represented by a polynomial curve-fitting 
technique as : 


Lua = Lug = 0.04765 * 18,—0.3129* If, +1.18« 13,—2.57* 12,+2.636* Im +0.15 


Figure (4-a.b) shows the Matlab and practical results for the starting-up condition of the 
generator terminal voltage. Figure (5-a.b) shows the Matlab and practical results of the steady- 
state phase voltage of the PMSG. From these figures, there is a very close agreement between 
the simulation and practical results, which indicates the validity of the proposed dynamic 
simulation model. Figures (6-7) show the simulation results of PMSG electromagnetic 
generated torque and output active power per phase of the generator. There is very close 
agreement between the theoretical and practical results due to all machine parameters, 
especially iron core loss, stray load loss, and the saturation effect are considered in the 
comprehensive dynamical model of the PMSG. 
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Figure 4-a,b Experimental result (a) and simulation result (b) of starting—up generator terminal 
voltage at the no-load condition. 
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Figure 5-a,b Experimental result (a) and simulation result (b) of the steady-state generator phase 
voltage 





Figure 6 Simulation result of electromagnetic torque Figure 7 Simulation result of active power per phase 


4. CONCLUSION 


A comprehensive dynamic model of PMSG 1s considered in detail in this paper taking the effect 
of iron core losses and stray load losses of the machine into account without increasing the 
number of differential equations of the model. The effect of iron core loss is represented by a 
resistance reflected in the stator circuit as a voltage drop. The stray load loss is represented by 
a series resistance connected in the stator circuit. The equivalent circuit of PMSG is modified 
to deal with all machine parameters without losing the accuracy of the machine performance 
determination. A series capacitance is connected in each generator terminal to obtain maximum 
output power from the generator when operates at synchronous speed and resonance conditions. 
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